Dyson-Schwinger equations with a parameterized metric 
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We construct and solve the Dyson-Schwinger equation (DSE) of quark propagator with a pa- 
rameterized metric, which connects the Euclidean metric with the Minkowskian one. We show, in 
some models, the Minkowskian vacuum is different from the Euclidean vacuum. The usual ana- 
lytic continuation of Green function does not make sense in these cases. While with the algorithm 
we proposed and the quark-gluon vertex ansatz which preserves the Ward-Takahashi identity, the 
vacuum keeps being unchanged in the evolution of the metric. In this case, analytic continuation 
becomes meaningful and can be fully carried out. 
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I. INTRODUCTION 

The true metric of space-time is certainly the 
Minkowskian other than the Euclidean. However, with 
the Minkowskian metric, in many circumstances, it is 
very difficult to perform calculations of physical quan- 
tities of our interests. Such difficulties are mainly 
caused by the kinematics which are only included in the 
Minkowskian system but not in the Euclidean one. Nev- 
ertheless, in the framework of path integrals in quantum 
field theory, the vacuums implied by the Minkowskian 
metric (with an infinitesimal e-term) and the Euclidean 
metric should be the same (see, for example, Ref. [1| 
or section II of this paper). Therefore, if we only fo- 
cus on some static quantities such as density, conden- 
sate and space- like Green function, Euclidean metric and 
Minkowskian metric are equivalent, giving the same ex- 
pectation value. However, if we also have interests in the 
physical quantities which involve in transport processes 
such as conductance and viscosity, the Minkowskian met- 
ric would be inevitable. In recent years, viscosity of 
strongly coupled quark-gluon plasma (sQGP) [2[ and the 
problem of energy losing during quark jet quenching pro- 
cess in sQGP have attracted a great deal of attention. 
These problems all involve transport processes and re- 
quire Minkowskian metric to be implemented [3| . 

How could we realize Minkowskian metric from an Eu- 
clidean system? The most popular way is the analytic 
continuation (Wick rotation) if we know the concrete an- 
alytic form of the Green function. We can easily ap- 
ply such an approach in perturbative calculations. How- 
ever, in many non-perturbative cases, for examples, the 
lattice QCD calculations [Jl and the Dyson-Schwinger 
equation (DSE) approach |5[ which are constructed in 
Euclidean space, we can hardly perform an analytic con- 
tinuation directly to the numerical results calculated in 



Euclidean space. Of course we might use these results to 
fit an analytic form in some speculations [g, 0] • For in- 
stance, a useful priori decomposition rules is the Kallen- 
Lehmann spectral representation. However, the problem 
arises from the lack of a systematic step-by-step program 
to efficiently assess and improve the accuracy of these 
speculations. 

In many models built in Euclidean space, another even 
more serious problem bothering us is whether the ana- 
lytic continuation of Euclidean Green function is mean- 
ingful and self-consistent with the principle of quantum 
field theory, namely whether the Minkowskian path inte- 
gral formalism and the Euclidean one lead to the same 
vacuum. 

In this paper, we will develop a new scheme, in the 
framework of DSE, which is free of all the problems men- 
tioned above. The paper is organized as follows. At first, 
in section II, we introduce the parameterized metric in 
the context of Quantum Field Theory (QFT) and show 
how it works as an effective parameter for the task of an- 
alytic continuation. In section III, the so-called vacuum 
problem, which bothers us in some models, is described 
explicitly. Then, we construct Dyson-Schwinger equa- 
tion for quark with the parameterized metric in section 
IV. The numerical algorithm and calculated results are 
given in section V. Finally, in section VI, we make a sum- 
mary and give some remarks. 



II. PARAMETERIZED METRIC AND ITS 
VALIDITY IN SOME SIMPLE CASES 



To introduce a parameterized metric, we should deal 
with the metric as a variable everywhere in QFT. For 
example, we consider the simple </> 4 theory. Its action 
can be expressed as 
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The covariant Minkowskian/Euclidean metric adopted 
in this paper is gff v — diag(— 1,+1,+1,+1) and g® v = 
diag{-\-l, +1, +1, +1) and the contravariant metric sat- 
isfies g^ x g\ v = 8%. Considering the requirements of 

»/— detg® = — i and \/—detgM, = 1, we could parame- 
terize the covariant metric as 
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with e G [0,7r). Then we have g(e — 0) = g E and 
g(e — > 7r) — > g M . Furthermore, one could easily prove 
that iS M [4>) = iS[g(n),4>] and -S E [4>] = iS[g(0),4>], so 
that the uniform path-integral formalism reads 



m 



JS[g(e),<t>] 



(3) 



Up to now, we have only found a simple way (namely, 
varying the parameter e) to connect gM with g E - In 
fact, there are infinite ways to realize this connection. 
Why ours is preferred, or what does it mean? To answer 
these questions, as an simple example, we construct the 
quantum theory of the 4 model from the Lagrangian 
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We will use the Hamiltonian formulation to quantize it. 
The momentum density conjugate to 4>{x) is 



BT 

<*) = TT = -g°°(e)\/-detg^(e) (0). 



(4) 



By using the relation detg M „ = goo = -mr, we obtain the 
Hamiltonian 

H= j d 3 x[ir<j ) -L}=e^')/ 2 I d 3 x[^7r 2 + ^m 2 ^+^ 4 }. 

'(5) 
At t = 0, we suppose that 4>(x) and n(x) are Hermite 
operators and satisfy the relation 



(x),Tr(y)} = i5 i3 \x-y), 

(x),0(y)] = [7r(f),7r(y)]=O : 



(6) 



and the equation of motion for Heisenberg operator is 
0(t) = J7(t)- x O(0)I7(t), (7) 

where the evolution matrix is 

U(t) = e- lHt = e - iHMi . (8) 

Here, Hm is the Hamiltonian operator in Minkowskian 
space as the same as that commonly used, and t = 
e l ( e ~^l 2 t. We should note that U(t) is not a unitary 
matrix since the H is not Hermite. So the operators 



cf>(t) and n(t) need not to be Hermite when t / 0, but 
it does not affect Eq. ([5]) to be held at t j^ 0. In fact, if 
(f)(t = 0) is Hermite, from Eqs. (g]), J7J) and © we can 
prove that 7r(i = 0) is also Hermite (note that (f>(0) is not 
Hermite). Moreover, one could easily prove that the cor- 
relation function G(x,y) in vacuum is time(i)-translation 
invariant. From these observations, we conclude that our 
construction is self-consistent. 

In addition, there are some other properties that we 
need for our purpose in this paper. Firstly, the structure 
of eigenstates of the H is independent of e. Secondly, the 
effects of the parameterized metric we adopt is equivalent 
to the Wick rotation to time axis in complex plane but 
nothing else. In fact, one prefers to regard Eq. (JHJ as the 
standard definition of Wick rotation and take it to real- 
ize analytic continuation of Green functions. However, 
in path-integral formalism, our method is more conve- 
nient. As we know that the path-integral formalism can 
automatically perform the time order operator, but the 
in/out states are unfixed in it. Nevertheless, the parame- 
terized metric undertakes an extra responsibility, besides 
the Wick rotation, that guarantees the in/out states to 
be those in vacuum (the lowest energy state of Hm)- It 
is true in our method with varying e G [0,7r). To see it 
more explicitly, we just need to show that 

!7(00)|a) OC e -HMcas(e/2)oo | Q > ^ ^ ^ 

It should be emphasized that Eq. © is definitely true 
except for the case with e — it, so that we could only 
take the limit e — > it to realize Minkowskian space whose 
effects are equivalent to the £-term trigger adopted com- 
monly (see, for example, Ref. [l|). Until now, we could 
say that this parameterized metric, used in path-integral 
formalism for analytic continuation of Green functions, 
is legitimate. (In addition, Eq. (J9j) means that the vac- 
uum state implied by g(e) is identical and independent 
of e. However, as we will see in the next section, such a 
statement is not always satisfied by models. We refer it 
as vacuum problem.) 

As a simple test of our method, we implement it to 
calculate the one-loop electron self-energy At first, we 
need conventions to simplify our representation. In this 
paper, real space-time variables are always thought to be 
the components in a contravariant vector x^ . Naturally, 
the real momentum variables are always thought to be 
those in a covariant vector p^ . The Fourier transform is 

/(*") = / ^y/-dotg^(e)e- i ^"f{p li ) , (10) 

/(p M ) = J d 4 x^-detg^( £ y p ^f(xn ■ (11) 

Since we are considering Fermion, the Dirac matrix 7^ 
should also be parameterized to satisfy 



{y(£),7^)} = 2<r(e), 
{l^),lu{e)} =2g liV (e), 
{7 M (e),7„(e)} = 2ff£(e) = 2#. 



(12) 



With the parameterized metric, the electron self-energy 
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fc 2 = g a 0k a k? . 

Using the Pauli-Villars regularization procedure, we have 
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where the functions A and B can be written explicitly as 



A(P 2 ) 



B(P 2 ) 
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which can be calculated by proper numerical integral. 

Note that £2 (p) can also be obtained through standard 
perturbative procedure, which reads 

a f 1 
s 2(p) = --r- / dx[2i(l - x)p + 4m e ] 

4tt7 

(l-x)A 2 



x In 



x(l — x)p 2 + xm\ 



(20) 



From Eq. (|2U)1 . we can also extract the structure func- 
tions A(p 2 ) and B(p 2 ), where p 2 should be understood 
as complex number but not that in Eq. (|14p . 

To illustrate the validity of the parameterized met- 
ric, we take the analytical structure of the mass function 
M (p 2 ) = B(p 2 )/A(p 2 ) as an example. The calculated re- 
sults in our presently proposed parameterized metric and 
that in the standard perturbative procedure are displayed 
in FigfTJ Looking over Fig[T] one can recognize defi- 
nitely that our presently proposed parameterized met- 
ric method is equivalent to the traditional perturbative 
calculations. 
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FIG. 1: (color online) Calculated result of the analytical 
structure of the mass function in our parameterized metric 
and the comparison with those given in the standard pertur- 
bative procedure, where p 2 — pe~ le . The solid line is that 
given by perturbative theory in the exact Minkowskian space 
(with e = tv in the expression of p 2 ), the dot line is also that 
calculated in perturbative theory but with e = 0.957T, the 
dash line is that obtained with our parameterized metric at 

£ = 0.957T. 



It should be noted that, in this example, the equiva- 
lence between our newly proposed method and the stan- 
dard perturbative calculation is somewhat trivial and our 
method is not economic. However, in non-perturbative 
frameworks, traditional analytic continuation does not 
work anymore (it requires us knowing the strict analyti- 
cal expression of the solutions) , instead, our method will 
become valuable and efficient. Anyway, such an equiv- 
alence indicates that presently proposed parametrized 
metric method is absolutely valid. Then, as we will see 
below, we can successfully apply this method to Dyson- 
Schwinger equation. However, before we do that, there is 
another problem deserving our attention, namely the vac- 
uum problem, because it will influence the proper trun- 
cation of Dyson-Schwinger equation. 



III. VACUUM PROBLEM 

To demonstrate the vacuum problem clearly, we first 
consider a free quark system with finite mass and chem- 
ical potential. In this case, the analytical expression for 
quark propagator is exactly known. The total influence 
of chemical potential /j, on the quark propagator is just 
a complex continuation of p (p — > pg — ifj,). Actu- 
ally, in many articles, /z is used as a single parameter 
for analytic continuation and not referred to as chemical 
potential. When /i > to, Fermi-surface will emerge from 
the origin of momentum space, and the vacuum will be 
changed. However, since there is no interaction, the in- 
fluence of the chemical potential /i>mon the vacuum 
can not alter the quark propagator's form. Now, we turn 
on some interaction which needs not to be strong (see, for 
example, Refs. 0-@). In this case, the quark propagator 
will include self-energy part which is dependent on the 
vacuum. Therefore, as the vacuum has been altered, the 
quark propagator with fj, > to is no longer an analytically 
continued Green function with po — > po — i\i. Up to now, 
we have not questioned the meaning of analytic continua- 
tion. We could just say that, beyond the mass threshold, 
chemical potential is no longer a valid parameter for the 
task of the analytic continuation [7] . 

However, in reality, there might be even worse situation 
when we try to use the parameterized metric g^(e) to 
realize the analytic continuation. As we will see, with the 
increasing of the parameter e, the range of p 2 = p^p^g^ 11 
{jpn G R) will become larger and larger in the complex 
p 2 plane. Considering the Minkowskian limit e — > n, the 
maximal range of the complex p 2 would be in the lower 
half-plane except the time-like real axis. Once a pole en- 
ters the range of p 2 , the situation becomes very similar 
to what happened in the case of large chemical potential 
fi > to. The vacuum will be changed from then on. Be- 
cause of the interaction, Green function is no longer an 
analytic function of p 2 . In this case, how could we de- 
termine the physical one naively? Can we still say that 
we have chosen a wrong way of analytic continuation as 
the same as that for chemical potential? In the sense of 
analytic continuation, we are indeed wrong. However, at 
the same time, we also see that the results of analytic 
continuation (if exists) is no longer corresponding to the 
true vacuum structure implied by the Minkowskian met- 
ric. It will result in a serious contradiction. On one hand, 
choosing the strict analytic continuation means we have 
to abandon the Minkowskian metric which is our original 
and real purpose! On the other hand, choosing a com- 
plete Minkowskian framework, we will lose the analyticity 
of Green function! Obviously, the root of this problem is 
the presence of singularities outside of the time-like real 
axis of p 2 plane. Actually, in this case, the Minkowskian 
vacuum is different from the Euclidean one. For exam- 
ple, if the quark propagator has a pair of complex poles 



with Im(z) ^ (see, for example, Ref. |To|). 

S(p) = - + , 

l l ' P + z l l ' P + z 

the chiral quark condensate reads 

d 3 p 



(21) 



< qq >= 2N C Z 4 



(2nY 



M\P\), (22) 



with the 3-momentum distribution function 



1 f°° 

h(\ P I) = j- / dp 4 tr D S(p) . 



(23) 



For our purpose, we just need to evaluate the Euclidean 
fi(0) and the Minkowskian / 2 M (0) 



fi(0) = C d Pi 



p\ + z 2 p 2 + (z*) 2 



(24) 



ff(0) = C idp 4 [ 



l pl~ z 2 + M+ p\- (z*) 2 + iO +i ' 

(25) 
After some direct calculations, we obtain the difference 
between the distribution functions of the condensate in 
the two metrics as 



/ 2 B (0)-/ 2 M (0) = 2^C^0. 



(26) 



Then the chiral condensate takes definitely different val- 
ues in the two metrics. 

Since qq is a gauge invariant local operator, different 
condensates (or its distribution) must come from different 
physical vacuums. It is easy for us to accept a different 
vacuum induced by the chemical potential /i > to. How- 
ever, in the framework of path integrals in quantum field 
theory, the vacuums implied by the Minkowskian met- 
ric (e — > 7r) and the Euclidean metric must be the same 
as we have already seen in section II. So, if a model is 
well defined, the contradiction mentioned above should 
not exist. It means that the singularities of the propaga- 
tor can not exist outside of the time-like real axis of p 2 . 
With this restriction, the parameterized metric can help 
us to perform the analytic continuation successfully and 
will indeed preserve the vacuum until g^(s) approaches 
infinitely the Minkowskian metric. 

In the following, we will try to establish such an ap- 
proach in the framework of the Dyson- Schwinger equa- 
tions 0. 



IV. QUARK'S DSE MODELS WITH THE 
PARAMETERIZED METRIC 



In this paper, the quantity of our interests is the quark 
propagator in momentum space 



Sip) = / d*Xy/-tetg» v J*- x ip.\Tq(x)q{f))\Sl) . (27) 



It satisfies the quark Dyson- Schwinger equation [5] 

S- 1 (p) = (v r -P + m) + E(p), (28) 

with 

S(p) = J V-detg^g 2 D^(k)^S(q)T:( qiP ) . 

(29) 
And according to its Lorentz structure, it is usually de- 
composed as 



or 



S- l (p) = i 1 .pA{p 2 ) + B(p 2 ), 



S{p) = i-f.pa A {p 2 ) + cr B (p 2 ) . 



(30) 



(31) 



Since we will make use of a model that is free of ul- 
traviolet divergence, here, we need not to discuss the 
issues related to the regularization and renormalization 
of Dyson-Schwinger equation [111 . Il2j . Anyway, if one 
adopts a model that has ultraviolet divergence, Pauli- 
Villars regularization would be a reasonable choice just 
as we have done in section II. 

Because the quark DSE involves effective (full) gluon 
propagator and quark-gluon interaction vertex which can 
not be determined by quark DSE itself, we need to choose 
some models for them. In Landau gauge, we have 



D^{k) = D{k)[g^{e) 



k^k v 



(32) 



with k^ = g^ v k v . And for the vertex r" Abelian as- 
sumption gives 



T a f0 (q,p) = ^-r il ( q ,p). 



(33) 



Finally, we need to choose some proper models for D(k) 
and T M (q, p) . In this paper, to preserve the vacuum in the 
evolution of the metric from Euclidean to Minkowskian, 
we have to adopt a self-consistent gluon propagator 
whose singularities do not appear outside of the time- 
like k 2 < axis. As a simple choice, we take 



g 2 D(k) = 4tt 2 D 



X 



(fc 2 +A) 2 



(34) 



in this paper. Since gluon is confined in reality, one 
might have instinctive suspicion on our choice in sense of 
confinement. Nevertheless, the lack of Kallen-Lehmann 
spectral representation is not the only way to exclude 
a particle from the physical Hillbert space. For exam- 
ple, in QED, the longitudinal photon indeed have legit- 
imate Kallen-Lehmann spectral representation, but it is 
still located outside of the physical Hillbert space. Of 
course, we are far away from the solution of QCD con- 
finement [13l4l5j | . People need to make much more efforts 
for this problem. More or less, it has gone beyond the 



main scope of our current discussions. Anyway, we will 
come back to this issue later. 

For the quark-gluon interaction vertex, the simplest 
choice is the so called rainbow approximation 



r^ = 7 ^ 



(35) 



In many practical calculations, the naive approximation 
has got triumphs. However, given the Abelian assump- 
tion, the bare quark-gluon interaction vertex violates 
the Ward-Takahashi identity (WTI) which is the con- 
sequence of the conservation of quark current. During a 
single-quark process, current's conservation is crucial for 
the unitary. Therefore, rainbow approximation could not 
promise that singularities of quark propagator can only 
appear at the time-like real axis (p 2 < 0). 

Beyond the bare (or rainbow) approximation, people 
have developed other vertex forms that satisfy the WTI. 
It has been shown that the so-called Ball-Chiu (BC) 
ansatz for the vertex 11611 is a successful form of them 
(see, for example, Ref. [17J)- It is expressed as 






with 



^Al f ,+ 1 ' {P +q) A A (p+q), l+i A B (p+q) fi , (36) 



A^ 
A B 



A(p 2 )+A(q 2 ) 

2 
A( P 2 ) - A( P 2 ) 

9 2 

£T — q A 

B(p 2 ) - B(p 2 ) 
p 2 — q 2 



(37) 
(38) 
(39) 



Since BC vertex fully satisfies WTI, we expect that sin- 
gularities of quark propagator could only appear at the 
time-like real axis (p 2 < 0) (similar conclusion has al- 
ready been drawn in Ref. 0). Then, we will take BC 
vertex in our practical calculation. 



NUMERICAL CALCULATION AND 
RESULTS 



A. Algorithm and its validity 

For numerical calculation, we extract at first the equa- 
tions for the functions A and B from the DSE of quark. 
After some derivation, we can write the equations about 
the A and B in the parametrized metric as 



A(p 2 ) = l+/e ^dq^dq 



|q| e A [A(p 2 ),A(q 2 ),p 2 ,q 



4M q 2 A 2 (q 2 ) + B 2 (q 2 ) ' 

(40) 



B(p») - /e-^ W lq| e ^ 2 )'^ 2 )'^ 2 ] 



4| go | q 2 A 2 (q 2 )+B 2 (q 2 ) 



(41) 



with 



9.4 = A(q 2 )V A F\ + A(q 2 )A A F; 
e B = B(q 2 )Z A F B + B(q 2 )A A F t 



B(q 2 )A B F A , 
-A(q 2 )A B F B , 



where q 2 — q^e %e + q 2 , and F\ B (i = 1, 2, 3) are func- 
tions of p 2 and q 2 because the angular integrals have been 
done. These simplified equations are the starting point 
of our further numerical calculations. 

Note that a given metric g^ u (e) determines a specific 
region of p 2 on complex plane which can be shown in 
Fig. [2] Then we can parameterize p 2 as pe~ lB with 
p € [0,oo) and 6 [0, e]. Obviously, when e ap- 
proaches 7r, the region of p 2 plus its conjugation will 
cover the whole complex plane except the time-like axis. 
It should be mentioned that, if we take the translation 
of po ~> Pa — ip to perform the analytic continuation, 
the range of (po — ip) 2 + p 2 is bounded by a parabolic 
curve. Even if singularities only appear at time-like real 
axis, such parabolic curve will inevitably encounter the 
singularities as p become moderately large. In this case, 
the available analytic information is only the inner part 
bounded by the parabola which can not include any sin- 
gularity. However, with the present parameterized met- 
ric, we can easily overcome such a problem. 




FIG. 2: [color online] Schematics of the area and the contour 
of the complex momentum p 2 (the lower part) and p 2 * (the 
upper part) with the parameterized metric g M "(e). 



With above representations, we obtain 



J 



^■n-^f*^ a 



iin(e - 6') \ ® A [A(pe- ie ), A(p'e- W ),pe~ l \ p'e 



sinO' 



B(pe~ ie ) = (e~ 1 ^ 2 ) / p'dp' / d6' 
Jo Jo 



4sine 



(p'e- 10 ') A 2 (p'e- 10 ') + B 2 (p'e- ie ') 



sin(s - 6') \ <d B [A(pe- 10 ), A(p ' e~ 10 ), pe~ l \ p'e 



-iB'l 



sine' I (p'e- w ')A 2 (p'e- ie ') + B 2 (p'e- 10 ') 



r 



(42) 



In principle, the above equations can be solved directly. 
In practice, to simplify numerical calculations, we con- 
sider ultraviolet boundary conditions for the functions A 
and B and the Cauchy integral formulas [lOj which read 



A(C 3 )=A*(C 1 ), B(C 3 )=B*(C 1 ), 


(43) 


A(C 2 )=1, B(C 2 ) = 0, 


(44) 


Z7rj JC1+C2+C3 l z 


(45) 



where the paths C\, C2 and C3 are those shown in Fig. [2] 
and F stands for the functions A and B. Then, in the 
integral equations, the only unknown function is F(C\) 1 
or F(pe- le ) which can be fully solved. Therefore, the 
original two-dimensional problem is reduced to a one- 
dimensional one. Such a reduction is very significant 
because it makes the numerical calculation with large e 
(— ¥ 7r) possible. 

To check the validity of our method, we adopt firstly 
a relatively small e = 37r/4 to solve the equation (with 



D = 0.25GeV 2 , X = 1-OGeV 2 and A = O.lGeV 2 ). The 
solutions of the functions A and B are illustrated in 
Fig. El Making use of these solutions, we can generate 
the data at 9 = with the Cauchy integral formulas. 
We compare then the generated data with the solutions 
obtained by solving directly the equations in Euclidean 
space (solving the equations after taking e = 0). The 
comparison is displayed in Fig. [4] The figure manifests 
apparently that the two groups of curves are exactly co- 
incident. It indicates that the validity of our method 
is in very high precision and the vacuum really remains 
against the evolution of e from to 0.75-7T. 

To clearly illustrate the self-consistency, we have to 
show that the quark propagator has no pole in the cor- 
responding area. We display then the obtained quark 
propagator (actually the absolute value of the cr Bl de- 
noted as I 0b |, is sufficient for this checking) in the whole 
range of p 2 in Fig. [5J It is apparent that there is indeed 
no pole located in the range of p 2 . 

In principle, for a given e, as long as the Green func- 
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FIG. 3: [color online] Calculated variation behaviors of 
A(pe~ le ) and B(pe~ le ) with respect to p at e = 37r/4. 
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FIG. 4: [color online] The comparison between the data gen- 
erated by our algorithm and those obtained directly in the 
Euclidean space. 



tions have no singularity in the involved momentum re- 
gion, our results would be reliable. As we have men- 
tioned above, the singularities of the quark propagator 
could only appear at the time-like axis, then the analytic 
continuation on the whole p 2 plane except the time-like 
axis could be achieved by the limit e — Y 7r. 



B. Almost Minkowskian propagator 

In this subsection, we discuss the quark propagator in 
the case that is very close to the Minkowskian by per- 
forming the calculations with increasing the £ to a value 
highly close to 7r. It should be noted that, the true limit 
(e = 7r) is very hard to achieve (and not necessary) in 
practice. At our current level, the largest e that can 
yield reliable results is 0.957T. In addition, we consider 
the dependence of the quark propagator on the coupling 
strength D and the screening width A. 

After having solved the quark DSE with our algo- 
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FIG. 5: [color online] Calculated variation behavior of the 
scalar part of the quark propagator |cr s | in the whole mo- 
mentum area in the case of e — 37r/4 (upper panel) and that 
at two special values of 8 (lower panel). 



rithm, we first analyze the behaviors of the functions 
A{p 2 ) and B{p 2 ) in the quark propagator against the pa- 
rameter e and with fixed parameter set D = 1.06GeV 2 , 
X = l.OGeV 2 and A = 0.5GeV 2 . The obtained results of 
the real part of the functions A and B at several values of 
e arc illustrated in the upper, middle panel of Fig. [6] re- 
spectively, and both the real and imaginary parts of the A 
and B at the maximal value of e we reached in the lower 
panel of Fig. [5] We can notice apparently from the figure 
that both the functions A and B have a peak at some 
complex momentum. Moveover, the peaks grow sharper 
and sharper with the increasing of e. Although we have 
to stop at £ = 0.957T because of numerical difficulty (the 
peaks make the numerical contour integral inexact as the 
e approaches much more to n), we can predict that the 
peaks will finally evolve into singularities. Namely, the 
analytic functions A(p 2 ) and B(p 2 ) both have poles at p 2 
(on the time- like momentum axis). 

To show the limit behavior of the imaginary part of 
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FIG. 6: [color online] Calculated results of the real part of 
the functions A and B in the quark propagator (Re(A) in 
upper panel, Re(_B) in middle panel) with several values of 
parameter e and those of the real and the imaginary parts of 
the functions A and B at s = 0.95-7T (lower panel). 



A and B as e — > n, we need the results with e much 
more close to it. Up to now, we can not say any definite 
things about the branch-cut structure of the quark prop- 
agator yet. However, there are some indications that the 
power of the divergence of A and B is less than one and 
there would exist branch-cut on the time-like axis. For 
instance, we find that the absolute values of A and B in- 
volve oscillation when p 2 < p 2 which seems to come from 
logarithmic divergence. Another evidence is the sym- 



metric behavior shown in Fig. [5] Since Fig. [5] is drawn 
in lnp 2 , the symmetric feature around the pole also sup- 
ports logarithmic divergence. However, this issue needs 
further investigation. 

Then, from the behaviors of the functions A and B, 
we can obtain the quark propagator. The obtained re- 
sults of the \a B | in the whole momentum region with the 
reachable maximal e (= 0.957r) is shown in Fig. [7] One 
could easily find from the figure that the Minkowskian 
quark propagator has one pole at p 2 , one zero-point at 
p 2 respectively. Noting that both the A and B are 
approximately constant at p 2 > p 2 which means that 
the dynamical mass function M{p 2 ) = B(p 2 )/A(p 2 ) also 
roughly remains constant, one can then recognize that, as 
the time-like momentum p 2 ss — M 2 (0) > p 2 , the struc- 
ture l/[p 2 + M 2 (p 2 )] will generate a simple pole. Ob- 

ned 
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FIG. 7: [color online] Calculated result of the variation be- 
havior of the function |er s | in the whole momentum area in 
the case of e = 0.957T. 

It should be emphasized that the simple pole of the 
Minkowskian quark propagator will likely lead to an un- 
confincd quark. Although the result is explicitly against 
the property of confinement, we tend to think that it is 
caused by the specific model for the gluon propagator. 
Actually, the gluon model adopted in this paper is not a 
rigorously confined one, because the linear confinement 
potential requires A — > 0. Then, we wonder what will 
happen when we decrease the parameter A and main- 
tain the mass gap M(0). The obtained results at several 
values of the A are presented in Table |U It is evident 
that p 2 approaches p 2 with the decreasing of A. It is a 
positive signal for our assumption that A — > will lead 
to p 2 — > pi. However, we could not prove this statement 



yet, since numerical difficulty prevents us from decreas- 
ing the A to infinitesimal. In this article, we leave this 
question open. However, given such an assumption, we 
could expect that, as pi —¥ p 2 , the simple pole and the 
zero point of the quark propagator will combine into a 
lower-order pole (fractional power or logarithmic diver- 
gence) which finally confines quark. 



TABLE I: Calculated results of the parameter dependence of 
the Minkowskian quark propagator 



D (GeV 2 ) 


A (GeV 2 ) 


vl (GeV 2 ) 


M(0) (GeV) 


Pi (GeV 2 ) 


1.06 


0.50 


-0.88 


0.11 


-0.012 


0.66 


0.30 


-0.58 


0.11 


-0.012 


0.25 


0.10 


-0.22 


0.11 


-0.012 



VI. SUMMARY AND REMARKS 

When we have interests with the quantities, such as 
conductance and viscosity, which involve transport pro- 
cesses, the Minkowskian metric would be inevitable. The 
traditional method to realize the Minkowskian metric is 
the so-called analytic continuation of Green function. In 
this paper, we argue that the meaningful analytic con- 
tinuation of Green functions requires an identical vac- 
uum. It should be true for any rigorous quantum field 
theory, but not always promised by models. Besides 
this problem, the widely used method for analytic con- 
tinuation in various models is based on the traditional 
K alien- Lehmann spectral representation. However, in 
some circumstances, especially for quark fields, the situa- 
tion would become more subtle and complicated. In vac- 
uum, since quark is confined, quark propagator may not 
have standard K alien- Lehmann spectral representation. 
More or less, such a difficulty may even exist in hot/dense 
quark matter where the quark degree of freedom becomes 
important. Therefore, it will invalidate the method of an- 
alytic continuation based on the K alien- Lehmann spec- 
tral representation. It has also been emphasized that, if 
we implement the translation po —* Pa ~ *M to perform 
the analytic continuation, the range of (p — i^) 2 + P 2 
is bounded by a parabolic curve. Even if singularities 
only appear at time-like real axis, such parabolic curve 
will inevitably encounter the singularities as /i become 
moderately large. In this case, the available analytic in- 
formation is only that in the inner part of the parabola 
whose top approaches the first singularity. 

In this paper, we have developed a new scheme which 
is free of all the problems mentioned above. A well- 
behavior effective gluon propagator and the Ball-Chiu 
vertex ansatz can promise that the singularities can only 
appear at time- like axis, thus promise an uniform vacuum 



under Euclidean and Minkowskian metrics. Furthermore, 
by establishing a parameterized metric g tJlfJ '(e) which con- 
nects the Euclidean metric with the Minkowskian one, we 
can directly solve the Green function in the full p 2 plane 
except the time- like axis without implementing any priori 
decomposition rules such as the K alien- Lehmann spectral 
representation. After confirming the correctness of our 
parameterized metric scheme in the simple (j) 4 model and 
in the calculation of electron self-energy at one-loop level 
and checking the validity of the algorithms for solving 
the quark DSE in the case of a relative small parameter e 
(= 0.757r), we continuously increase the e and obtain the 
almost Minkowskian quark propagator with e — 0.957T. 
We find that the Minkowskian quark propagator has one 
simple pole and one zero-point at pi and pi respectively. 
By observation on the dependence of the quark prop- 
agator on the screening width A in the effective gluon 
propagator, we could expect that, as A — > 0, the simple 
pole and the zero point of quark propagator will combine 
into a lower-order pole (fractional power or logarithmic 
divergence) which finally confines quark. However, the 
exact structure (for instance, branch-cut) has not yet 
been identified due to the limited e. In hot/dense strong 
interaction matter where the quark degree of freedom be- 
comes more significant and practical @, 0, 0, H, HH, [l9j , 
our method need further modifications which is under 
way. In this case, an interesting question is whether 
the so-called deconfined quark has the standard Kallen- 
Lehmann spectral representation. Anyway, an opposite 
answer can only come from the non-perturbative nature. 

In addition, if one makes use of a gluon propagator 
model, whose ultraviolet behavior being consistent with 
the perturbative result, the loop integral in quark DSE 
will be divergent. In this case, proper regularization 
would be inevitable. To be honest, we have no idea about 
how to take dimensional regularization in the numerical 
calculation at present stage. Since p 2 is a complex num- 
ber in the calculation, the regularization set by p 2 < A 2 
is also invalid. One might take \p 2 \ < A 2 alternatively. 
However, it is a non-analytic cut-off and inconsistent with 
our method. In our view, Pauli-Villars regularization 
procedure could be a reasonable choice, as we have taken 
in section II, which can overcome the difficulties men- 
tioned above. The related investigation is in progress. 
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